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\^ Abstract 

(N 

^ In sneutrino liybrid inflation tlie superpartner of one of the right-handed neu- 

trinos involved in the seesaw mechanism plays the role of the inflaton field. It 
obtains its large mass after the "waterfall" phase transition which ends hybrid 
inflation. After this phase transition the oscillations of the sneutrino inflaton field 
'""' may dominate the universe and efficiently produce the baryon asymmetry of the 

^^ universe via nonthermal leptogenesis. We investigate the conditions under which 

OO inflation, with primordial perturbations in accordance with the latest WMAP re- 

O 

[^ suits, as well as successful nonthermal leptogenesis can be realized simultaneously 






q 

o 
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within the sneutrino hybrid inflation scenario. We point out which requirements 
^^ successful inflation and leptogenesis impose on the seesaw parameters, i.e. on the 

Yukawa couplings and the mass of the right-handed (s) neutrino, and derive the 
^ predictions for the CMB observables in terms of the right-handed (s)neutrino 

S^ mass and the other relevant model parameters. 



1 Introduction 

The paradigm of cosmic inflation [H [21 E] (for recent reviews see e.g. [H [S]) has proven 
very successful in resolving the flatness and horizon problems of the early universe and in 
explaining the absence of relics from early phase transitions. However, the connection 
to particle physics is still unclear. One possibility to establish such a connection is 
provided by sneutrino hybrid inflation [12j, where the superpartner of one of the right- 
handed neutrinos involved in the seesaw mechanism [7] plays the role of the inflaton 
field. In sneutrino hybrid inflation, a large vacuum energy density is present which 
drives inflation and the sneutrino direction in field space has an almost flat potential 
suitable for slow-roll inflation. The right-handed (s)neutrinos obtain their large masses 
after the "waterfall" phase transition which ends hybrid inflation. Inflation in this 
scenario is closely linked to the physics generating the small neutrino masses via the 
seesaw mechanism. 

Another attractive connection between the seesaw mechanism and early universe 
cosmology is the possibility of generating the observed baryon asymmetry of the uni- 
verse via the out-of-equilibrium decays of the right-handed (s) neutrinos in leptogene- 
sis [6] (for recent reviews see [HI |9j). For calculating the produced baryon asymmetry, 
the knowledge of the phase of (p)reheating after inflation is in general mandatory, 
since it may lead to the nonthermal production of right-handed (s) neutrinos and since 
it determines the reheat temperature, which in turn governs the possibility of thermal 
(s)neutrino production. In most inflation models the nonthermal (s)neutrino produc- 
tion must arise from the decays of the inflaton field. In sneutrino hybrid inflation, 
on the other hand, the inflaton itself is a right-handed sneutrino, which means that 
this intermediate step is skipped and the sneutrino inflaton field after inflation may di- 
rectly dominate the universe and, when it decays, most efficiently produce the baryon 
asymmetry and reheat the universe. 

In previous works, leptogenesis after sneutrino inflation has been studied in the 
context of chaotic sneutrino inflation fiUl [TT] which however requires a quite heavy 
sneutrino with a mass of about 10^^ GeV and correspondingly very small Yukawa 
couplings in order to realize a low reheat temperature as suggested by gravitino and 
similar constraints in supersymmetric cosmology. Furthermore, chaotic sneutrino in- 
flation with a quadratic potential for the inflaton gives rise to a comparatively large 
tensor-to-scalar ratio of r ~ 0.16. On the other hand, sneutrino hybrid inflation [T2], as 



typical for hybrid-type inflation scenarios fi3\ [HI [151 [IB] , predicts a much smaller ratio 
r < 0.01 and is thus clearly distinguishable from chaotic sneutrino inflation by future 
observations (e.g. by the Planck satellite). Recently, it has been shown that sneutrino 
hybrid inflation [12] belongs to a wider class of hybrid-like inflation models, dubbed 
"tribrid inflation" in [T7] , which are very suitable for being embedded into supergravity 
(SUGRA) theories with the SUGRA ?7-problem solved by either a shift symmetry r\ [2T] 
or a Heisenberg symmetry |22] in the Kahler potential. While the sneutrino was a 
gauge singlet in [12], it has been demonstrated in [23] that it may be embedded into a 
Grand Unified Theory (GUT) representation, e.g. into a 16-plet of SO(IO), establish- 
ing a possible link between sneutrino hybrid inflation and left-right symmetric GUTs. 
Nonthermal leptogenesis after sneutrino hybrid inflation, on the other hand, was only 
briefly discussed in [12] for an example set of model parameters. 

In this paper, we therefore investigate in detail the conditions under which inflation, 
with primordial perturbations in accordance with the latest WMAP results, as well as 
successful nonthermal leptogenesis can be realized simultaneously within the sneutrino 
hybrid inflation scenario. We point out which requirements successful inflation and lep- 
togenesis impose on the seesaw parameters, i.e. on the Yukawa couplings and the mass 
of the right-handed (s) neutrino, and derive the predictions for the CMB observables in 
terms of the right-handed (s) neutrino mass and the other relevant model parameters. 
Our results are meant as a guideline for the construction of explicit particle physics 
models incorporating sneutrino hybrid inflation and baryogenesis via nonthermal lep- 
togenesis. 

The paper is organized as follows: In section |2] we introduce the sneutrino hybrid 
inflation scenario in a simple setup. Section [3] is dedicated to the inflationary phase and 
the predictions for the CMB observables. In section |4] we discusses the reheating of the 
universe after inflation and the production of the baryon asymmetry of the universe 
via nonthermal leptogenesis. We conclude in section |5] by combining our results from 
inflation and leptogenesis to highlight the preferred ranges of the model parameters. 

2 Framework 

We will discuss sneutrino hybrid inflation and subsequent baryogenesis via nonthermal 
leptogenesis in an extension of the minimal supersymmetric standard model (MSSM) 



^In the context of chaotic inflation, shift symmetry has been used e.g. in [rSllTOlEO] . 



with conserved R-parity, where three additional right-handed (s) neutrinos acquire large 
masses after the "waterfall" phase transition ending inflation. The superpotential defin- 
ing our framework is given by 

W = WussM + iVuh N' k e^'Li + ^ {N'^fH^ + nS{H^- Kp) + . . . , (1) 

where the iV* (the index i = 1,2,3 denotes the different generations) are gauge sin- 
glet superfields describing the heavy right-handed (s)neutrinos and where the reduced 
Planck scale is given by Mp ~ 2.4 ■ 10^*^ GeV. The canonically normalized imaginary 
partqj A^* of the respective scalar components are infiaton candidates as will be de- 
scribed below. L and h are SU(2)2,-doublet superfields which contain the standard 
model leptons and up-type Higgs, respectively. The Yukawa coupling term of N with 
the Higgs and lepton doublet, i.e. the second term in Eq. ([I]), allows to identify N with 
the right-handed neutrino superfield. 

H and S are two additional gauge singlet superfields. Here the canonically nor- 
malized real part H of the scalar component of H is the so-called "waterfall" field 
responsible for ending inflation. The F-term of 5*, the so-called "driving superfield", 
provides the large vacuum energy density that drives inflation. The scalar component 



of S is flxed at zero during inflation by SUGRA corrections (cf. section 3.3) and does 
not affect the inflationary dynamics. Furthermore, we assume Xa and k to be real 
coupling parameters for simplicity. 

The form of the superpotential Eq. ([I]) is motivated as follows: The latter two terms 
generate the scalar potential suitable for inflation. In the false vacuum with large values 
of A^' and H stabilized at zero, the large vacuum energy Vq = h?M^ drives the quasi- 
exponential growth of the scale factor in inflation. Once the slow-rolling flelds A^* fall 
below a critical value, the negative contribution to the squared mass of H from the 
term kS{H'^ — M^) starts dominating over the positive contribution from the terms 
-^ (iV*)^if^. Therefore, H becomes tachyonic which triggers the "waterfall" ending 
inflation as H acquires a non-zero vacuum expectation value (vev). 

After inflation, close to the global minimum of the potential where N^ k. Q and 
H ^ \/2M and where the large vacuum energy contribution vanishes, the flelds A^* 
and H perform damped oscillations accounting for a matter dominated universe. The 



^At this point we will assume that inflation proceeds along the imaginary direction of the complex 
scalar sneutrino field. We will see later that this can be obtained by a shift symmetry in the Kahler 
potential, which protects this direction against the SUGRA ry-problem as was demonstrated in [3T]. 
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Figure 1: Typical F-term scalar potential resulting from the model of Eq. (fTl) for the scalar components 
of S, L, h set to zero. For the plot, we have used example parameters k — X — 1 and AI = Mp. 



field which decays last and finally dominates the universe is generically the right-handed 
sneutrino with the smallest mass and smallest Yukawa couplings, since it only decays 
via the second term in Eq. (jTl) proportional to {yu)ij- This decay reheats the universe 
which thus enters its radiation dominated epoch. For illustration, we have plotted 
a typical scalar potential resulting from the scenario described above in Fig. [1} For 
further details see e.g. Ref. ^H\ . 

The first three terms in Eq. (jTl) describe the MSSM with masses for the additional 
right-handed neutrinos generated after infiation. In particular the term -^{N^YH"^ 
generates mass terms for the heavy (s)neutrinos as H acquires its non-zero vev. The 
vev of H in the true minimum is governed by the fourth term in Eq. ([I]). In a realistic 
scenario, we would expect inflation to end by a phase transition, i.e. the H field to be 
a non-singlet under some symmetry group [j A realization within the context of GUTs 
can be found in Ref. [23], and the idea of relating the "waterfall" of hybrid infiation 
to the breaking of a family symmetry was discussed in Ref. [21]. For simplicity, we 
keep H a gauge singlet here. Finally, the ellipsis represent possible higher dimensional 
operators. We note that a superpotential of the form given in Eq. ([I]) can be realized 
in an explicit model with discrete symmetries, as has been shown in [211 123]. These 
discrete symmetries distinguish between the gauge singlet fields of Eq. ([I]). 



•^In this case, the terms H^ might be replaced by a different combination of fields, i.e. by H1H2, 
in order to form a singlet under the symmetry group. The symmetry could, e.g., be a unifying gauge 
symmetry or a family symmetry. 



The parameters which appear in Eq. ([I]), and which govern our framework, can be 
understood as follows: 

• The phase transition scale M is the vev of the scalar component of the H 
superfield after the phase transition ending inflation and is the mass scale relevant 
for inflation. 

• The parameter An determines the seesaw scale which corresponds to the mass 
of the lightest right-handed neutrino 7x1^1 = 2j^M'^ in the true vacuum of the 
theory. 

• The vacuum energy parameter k fixes the the vacuum energy density Vq = k^M^ 
with regard to the phase transition scale. 



• The effective Yukawa coupling yi = y {yuvl) 11 is linked to the light neutrino 
masses. 

We will be working in a SUGRA framework with SUGRA corrections stabilizing the 
scalar components of S, L and h during inflation and with a symmetry in the Kahler 
potential that guarantees tree-level fiat A^* directions. This solves the ?7-problem of 
SUGRA inflation. The details of such a SUGRA framework are discussed in section [3731 
However, to illustrate the underlying physics more clearly we first focus on a global 
supersymmetry (SUSY) model and take the features mentioned above for granted. 

In order to produce the CP-violation necessary for leptogenesis we work with three 
(s) neutrino generations. Assuming that the right-handed neutrinos are strongly hierar- 
chical, i.e. one of them is significantly lighter than the other two, the scalar components 
of the latter superfields can be stabilized at their minima before the final 60 e-folds of 
inflation begin. Thus the time evolution of the lightest sneutrino controls the relevant 
slow-roll dynamics and it can therefore be identified as the inflaton. On the other hand, 
the outcome of leptogenesis is governed by the sneutrino with the smallest decay rate. 
This implies a comparatively small mass and small Yukawa couplings. In the following, 
we shall concentrate on the case where the lightest sneutrino drives both inflation and 
leptogenesis. Hence, the three generation model can be simplified to an effective one 
generation model in the right-handed neutrino sector, with the only remnant of the 
other two generations being a non-vanishing CP-asymmetry necessary for leptogenesis. 
We can thus concentrate on i = 1 in Eqs. ([I]) and we denote the relevant inflaton 
direction by A^ = A^^ and the respective coupling constant by A = An. 



3 Inflation 

Based on the framework described in the previous section, we now have a closer look 
at the inflationary dynamics in our model. Furthermore, we derive restrictions on the 
model parameters from the requirement of successful inflation and the latest observa- 
tional data. We start with a short introduction to slow-roll inflation and then discuss 
a realization of the model of section |2] in a globally supersymmetric context. We then 
refine this discussion by including SUGRA effects and close the section by listing the 
inflationary predictions from our model and comparing them to the latest observational 
data. 

3.1 Short Overview 

A common way to realize inflation is the so-called slow-roll paradigm, where a classical 
scalar field with a strong dominance of its potential energy over its kinetic energy 
V ^ £kin drives the accelerated expansion of the universe. At the same time, the 
quantum fluctuations of the inflaton field can account for the metric perturbations 
which give rise to the small scale CMB anisotropics. Inflation ends when the slow-roll 
conditions are violated, i.e. when the slow-roll parameters parameterizing the scalar 
potential and its derivatives 

e=^(yj, V = M],-, e = M^,^^, (2) 

become of order one. Here, a prime denotes derivative w.r.t. the inflaton field A^. 

In the slow-roll approximation, when e ^ 1, |?7| <^ 1 and ^^ -C 1, the equation of 
motion of a homogeneous (classical) scalar field 

N + 3nN + V' = 0, (3) 

simplifies to 

3HN = -V'. (4) 

Here T-C denotes the Hubble expansion parameter. 

Models of inflation typically predict the power spectra of the gauge invariant scalar 
and tensor perturbations at the time when the relevant fluctuations exited the horizon, 
roughly A/'e ~ 50 — 70 e-folds before the end of inflation. The amplitude of the scalar 
perturbations A^, the scalar spectral index n^, the running of the scalar spectral index 



Usi the tensor-to-scalar ratio r and the tensor spectral index nt can be estimated in 
terms of the potential and the slow-roll parameters [1] as 

1 1 V^ 



K^ 



(5) 



Ml, 127r2 {V'Y ' 
TT-s^l — 6e + 2r7, 
as ~ 16er/-24e^-2^^ 
r ~ 16e , 
nt ~ -2e, 

where these expressions have to be evaluated at the field value A^ = N{Mc) when the 
relevant scales leave the horizon. This value can be computed from Eq. ^. 

In order to test our model against observations, we compare the predictions from 
Eqs. rtsl) to the experimental data obtained from the 7 year WMAP survey combined 
with measurements of the baryon acoustic oscillations (BAO) [25] and measurements 
of the present value of the Hubble parameter Hq [26] using the six parameter ACDM 
fit [2Z|, which are given by 



0.951 < n, < 0.975 
0.939 <n,< 0.987 



(68 % CL) , 
(95 % CL) , 



(6) 



A^ = (2.441+°-°^«) ■ 10-^ 



-0.092; 



3.2 Realization in Global Supersymmetry 

In the model described in section |2| we assumed a symmetry in the Kahler potential 
guaranteeing a fiat N direction (imaginary direction of the scalar component of A^^) at 
tree- level (see also section 3.3[ ). Thus loop corrections must be taken into account and 
these can indeed generate a small slope as required for slow-roll infiation. According 
to [2B] , the Coleman- Weinberg one-loop effective potential is given by 



V. 



loop 



647r2 



STr 



A^Mln 






(7) 



with Ai denoting the mass matrix of the theory and Q a renormalization scale. The 
A^-dependent bosonic and fermionic mass terms generating a slope for the inflaton via 
the loop potential can be calculated from the scalar F-term potential 

2 



v.=j: 






!>_!.$ 



(8) 



and the fermionic mass matrix 






(9) 



$* denote the superfields of the theory, $* the respective scalar components. The 
relevant (i.e. A^-dependent) contributions to the loop potential are give by the H, U 
and h mass terms u 



^^K:r = k^kT = {rn}^:r = ^N^^liv^jiji , 



AS)^2 _ (^{P)\2 _ (^{F).2 _ ^ Ar2 V^ |/„, ^ |2 



(-SV = (-SV = K?r = ^AT^iMi.p , 






{w}^'Y = 2k'M'{x-1) 



(10) 



{m!'Pf = 2^M\, 



with 






'p 



Note that the U and h terms in the supertrace vanish since the degeneracy in the 
respective fermionic and bosonic masses leads to a cancellation of these contributions. 
Embedding this model in SUGRA provides the necessary stabilization of the scalar 
components of the U and h superfields during inflation and removes this degeneracy. 
However in the parameter range of interest, the contribution of the U and h terms to 



the loop potential turn out to be negligible (see also section 3.3). In the following, we 
fix the renormalization scale to Q = \/2kM^ which is the order of magnitude of the 
SUSY breaking scale. In our model, inflation ends when the if -field destabilizes at the 
critical value A^'^ characterized by m\j = and thus x = 1, 

(Ar^)2 = y2^MMp. (12) 

A 

We can now determine the observables describing the CMB fluctuations given by 
Eqs. (tsl), thus obtaining expressions depending on the phase transition scale M, the 

■^Here the index (S) (for scalar) denotes mass terms of the real parts of the complex spin-0 compo- 
nents of the superfields whereas the index (P) (for pseudoscalar) marks the mass terms corresponding 
to the purely imaginary parts. The index (F) marks the mass terms of the fermionic components of 
the superfields. 

9 



seesaw scale rriN ~ A and the vacuum energy parameter k. With x as defined above, a 
Taylor expansion in 1/x (with 1/x < 1 because A^ > N'^ during infiation) yields 

Vioop ^ ^-^ In a; . (13) 

Inserting this into the equation of motion Eq. (|4]), with Ji approximately constant, 
gives the value for N ai JV e-folds before the end of inflation: 

N\M) = {N^y + '-^ Mj. . (14) 

With this, the inflationary predictions of Eqs. (IS]) are given by 

A: 



3k2M|, ' 

+ 4^ J vr2Ar2 ' 
a ^-,l + ^ + ^^^!!M (15) 






K^Ml 



3.3 Embedding in Supergravity 

We next consider a possible embedding of our model in SUGRA. We focus on a Kahler 
potential with the ?7-problem [151 129] resolved by a shift symmetry [21j in the inflaton 
direction 

K=\s\'+\H\'+\h\^+Y.\ {^' + i^yf + E \^'\' +^\sm^+.... (16) 

i j P 

The Kahler potential can be seen as a general expansion in the superfields of the 
theory with the additional feature of a shift symmetry which guarantees tree-level 
flat directions for the imaginary parts of the scalar components of N^ and thus possible 
inflaton directions. The A^-dependent mass terms generating a slope for the inflaton via 
the loop potential can be calculated from the scalar F-term potential and the fermionic 
mass matrix as before by 

Vf = e^ \k'W,WDjW* - 311^1^1 

L J <j.^$ ^^^^ 

{MF)ij = e^'/^Wij + KijW + K,Wj + K,Wi + K.KjW - K^'K,.{DkW) 



<J>— 5><3> 
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The relevant contributions to the loop potential are given by the H , U and h mass 
terms which now obtain SUGRA corrections: 



[m 



(5)n2 



H 



(m 



iP)\2 
H ) 



2k'^M^ 



2k'^M^ 



X 



1 



X + 1 + 



M 
M 



(1 - ksh) /2 



t^sn) /2 



m 



2 K^M'^x , 

AP)\2 _ 



H ) 

(-5)^2 



M^ 



\ni 



[m 



My 2 



|2 



I N'- j:\iv. 



l^)lj\ 5 



;i8) 



m 



[m 



(5)n2 
(i^)N2 



[m 



J 

{P)■^2 
Li ' 



^^^+^iv^i(..).r 



N\y, 



y)i3\ 



Comparing these expressions to the mass terms calculated in section 3.2 we note 



some important points. A second mass scale, the scale of the SUGRA mass split- 
ting kM'^ jMp^ has appeared. However this scale is much smaller than the SUSY mass 
splitting scale -\/2kM and thus we shall keep the latter scale as the renormalization 
scale. The additional mass splitting implies that the h and U contributions no longer 
cancel. However since the mass splitting is small compared to the SUSY mass splitting 
appearing in the H mass terms and since the remaining parts of the h and U mass 
terms are proportional to |(yjy)ij| these additional contributions to the loop potential 
are negligible for y\ < 10~^. We will see later that this easily holds in our model. 
Furthermore, a new parameter has appeared in the loop potential: 

• The SUGRA correction parameter ksh controls the SUGRA corrections to 
the loop potential, ush = 1 recovers the phenomenology of global SUSY. 

3.4 Predictions 

As in the globally supersymmetric case, predictions for observables describing the CMB 
spectrum can now be obtained by solving Eq. ^ and evaluating Eqs. ^ at the time 
when the CMB fluctuations exited the horizon. In the SUGRA scenario, this was done 
numerically for ^4 = 60. Fixing the phase transition scale M by the experimental 
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value for the amplitude of the CMB fluctuations A^, the behavior of the spectral index 
Ug, its running a^ and the tensor-to-scalar ratio r is shown in Fig. |2] Interpreting the 
results visualized in Fig. |2] and enforcing the experimental bounds of Eqs. (|6]) implies 
restrictions on the model parameters. 

• The phase transition scale M is fixed to M ~ 0.0032 MpC:i%- W^ GeV with a 
slight deviation in the region of large SUGRA corrections. This is consistent with 



the global SUSY calculation (from Eqs. (14) and (15)) which gives M^ ~ SA^/AQ 
for N{M^) > A^^ 

The width of the band in Fig. |2] is given by the variation of the vacuum energy 
parameter k. A priori we would expect k to be an 0{1) parameter, thus we shall 
assume 0.5 < k < 2. In Fig. [2| larger values of n are associated with SUGRA 
corrections becoming relevant at smaller values of irtM- In particular the tensor- 
to-scalar ratio r is quite sensitive to n with k = 2 leading to comparatively larger 
r~ 0(10-2). 

Fig. |2] also demonstrates the effect of the SUGRA correction parameter ksh- 
The respective quantities are marked in black for ksh — 1 = which corresponds 
to the globally supersymmetric limit and in blue (red) for ush — 1 = —1 (+1) 
which corresponds to turning on the SUGRA corrections in the H mass terms 
with positive (negative) sign. In the considered SUGRA context the value of ush 
is a priori undetermined. Thus we would in general not expect to find global 
SUSY restored, which would correspond to ksh exactly equal to one. 

The second parameter controlling the effect of the SUGRA corrections is the see- 
saw scale mj\f = 2XM'^/Mp. Fig. ^ shows that these corrections are suppressed 
for small rriN, i-e. the observables are independent of ksh for small values of the 
seesaw scale and the model predicts (for M fixed by the experimental value of A^) 

0.98<n, <1, 3-10-^<a, <0, r < 0.013 . (19) 



Note that in this case Eqs. ( 15 ) hold. For very small values of the seesaw scale ttin, 
the spectral index Ug approaches 1, which is not preferred by the latest WMAP 
data. On the other hand, all solutions with ksh 7^ 1 leave the experimentally 
preferred region for the spectral index at large values of rriN. In combination, we 
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Figure 2: Loop potential and predictions for CMB observables for ngu — 1 = — 1, Oand 1 (blue, black 
and red), (a) Loop potential for k. = 0.5, M = 0.0032 Mp, m^ = 2.5 • lO" GeV. (b) - (d): Spectral 
index, running of spectral index and tensor-to-scalar ratio. The width of the bands is given by the 
variation of the vacuum energy parameter k = 0.5 — 2. On the left border of plots (b) and (c), k — 0.5 
corresponds to the upper set of lines, and in plot (d) to the lower set of lines. For labeling the x-axis, 
the phase transition scale was set to M = 0.0032 Mp. The 95% and 68% CL experimental bounds 
from Eqs. (l6| arc marked by dashed lines in (b). 
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find the preferred regions 2 ■ 10 GeV < rriN ^ 7 • 10 GeV for ksh — 1 = +1 and 
2 ■ 10^° GeV < m;v < 2 ■ 10^^ Q^y fo^ f^^^ - i = _i, respectively. Q 

Motivated by the above results, we take the phase transition scale to be fixed at M ~ 
8 • 10^^ GeV and concentrate on the parameter ranges 

2 ■ 10^° GeV < mjv < 7 • 10^2 QgY ^ 0.5 < k < 2 , \ksh - M > O-l , (20) 

in the further discussion. This yields 

- 0.0004 < a, < 0.0002 , r < 0.015 , (21) 

for the running of the spectral index and the tensor-to-scalar ratio. 

4 Reheating and Leptogenesis 

After the end of the inflationary epoch the homogeneous classical fields and their quan- 
tum fluctuations evolve according to their respective equations of motion. The universe 
enters a matter dominated regime until the decay of heavy particles and the thermaliza- 
tion of the light particles result in the total energy density being dominated by radiation 
(Fig. [3]). This out-of-equilibrium decay of heavy particles can furthermore produce the 
necessary lepton asymmetry. In the following we will study these processes for the clas- 
sical fields (reheating) and briefly comment on possible effects originating from their 
fluctuations (preheating). To this end we will start with the equations of motion for 
classical scalar fields, justify a simplification to Boltzmann equations and finally derive 
analytical expressions for the generated baryon asymmetry and the reheat tempera- 
ture. We finish by commenting on preheating via parametric resonance and tachyonic 
preheating in this context. 

4.1 Classical Field Dynamics after Inflation 

The equations of motion for the scalar fields can be obtained by adding a phenomeno- 
logical decay term [30j to Eq. dsl) thus giving 



+ 3?/0 + \/'(0) + r0 = O with (f) = {N,H}. (22) 

Equivalently, we obtain 5 • lO""* < A < 0.13(0.043) for ksh - 1 = +1 and ksh - 1 = -1, 



5 
respectively. 
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Adding a Boltzmann equation for the quickly thermalizin^ ultra-relativistic particles 
and the Friedmann equation, we arrive at a closed set of differential equations: 

iV + 3HiV+|^ + rjviV = 0, (23) 

■ dV ■ 

H + 3nH+-- + THH = 0, (24) 

oH 

PR + iTipR-rNPN -ThPh = 0, (25) 

{PN + PH + PR) = n^ , (26) 



3M|, 



with pr, ptv and pn denoting the energy densities of the ultra-relativistic particles, 
the A^-field and the i7-field respectively with Pn + Ph = N^f^ + H^/2 + V{N,H). 



Having solved Eq. (23)- (26), the lepton number density n^ can be calculated from the 
Boltzmann equation 

riL + 3'HnL = eiTN \-^3^h , (27) 

ruN^ ruNs 

with the CP-violation per (s)neutrino decay e for a hierarchical spectrum of right- 
handed neutrinos bounded by [331 EH [35] 



3 ^/Aml^^mN, .^„. 

e* < o 1-^ ■ I28j 

OTT (t>)^ 

Here {v) denotes the vacuum expectation value of the up-type Higgs. The lepton asym- 
metry is typically normalized to the entropy density s = 27r^(7*T^/45 [3] with the effec- 
tive number of degrees of freedom g^: = 915/4 for the MSSM particles. The asymmetry 
ul/s is transferred to the baryon sector via sphaleron processes ub = rfn^L where 
C is a number 0{1) depending on the field content of the model and the temperature 
Tsph when the sphalerons leave equilibrium. In the MSSM C ~ 1/3 [9J. The quantity 
measured today is ?7 = ^^ which can be calculated from the results above using the 
current conversion factor s = 7.04 n-y [3j. The second important physical quantity in 
the theory of reheating is the temperature of the universe when the universe becomes 
radiation dominated {Tn ~ H), the so-called reheat temperature. It can be calculated 
from the results above using T'^ = p^- 30/(5'* vr^). 



Eqs. (23) - (27) assume that both the A^-particles and the if -particles decay into 



ultra-relativistic particles with the respective decay rates T^ and Th- We will now 



^In principle, thermalization in the MSSM could be delayed if MSSM flat directions obtain large 
vevs (see e.g. [311 132]). However in our scenario, this is not the case since only the N direction is 
protected against large SUGRA corrections. 
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describe a possibility how to evaluate these quantities in our framework. At the be- 
ginning of the reheating phase 1-L ^ T^i ^h holds which implies that the decaying 
particles are damped predominantly by Hubble expansion, not by decays, and the pro- 
duced ultra-relativistic particles are strongly diluted. The decays become significant 
for t ~ 1-L~^ ~ min{r]^^,r^^}. At this stage it is safe to assume A^ <^ 1. In this limit 
the respective decay rates derived from Eq. ([I]) are 

-, ^ hm ^Mlk^^ = (Mill! 2^M- , (29) 



N^o A-K 2-n M. 



p 



Th ^ lim min{rH^^3jv3, ^m-^hL} = -^M^mm{{y^yl)33, ^l } , (30) 

with the sneutrino A^ decaying directly into lepton and Higgsino or slepton and Higgs 
and the H particles decaying predominantly into the heaviest fermionic neutrino (as- 
suming this is not strongly suppressed by kinematics) R which then in turn decays into 
lepton and Higgs or slepton and Higgsino. 



Note that the Boltzmann equations (25) and (27) imply a splitting of the total 



matter energy density pu into pn and p^, which is not straightforward if the respective 
degrees of freedom are highly coupled. However, since F^v ^ T h in our setting |j in 
the preferred region of parameter space (see section [s]), any radiation energy density 
produced by //-decays will be strongly diluted during the following matter dominated 
phase governed by oscillations of the sneutrino. With pM ~ Pn shortly after the 
end of inflation due to the strong damping of the i/-field (see below) we can thus 



substitute (25) and (27) by 



Pr + ^UPr-TnPm^Q, (31) 

rlL + 'iV.nL-eiTN^^, (32) 

without introducing a significant error for the finally produced radiation density. The 
effect of this approximation on the Hubble expansion rate is negligible since the Fried- 
mann equation is predominantly governed by pM for t < T~^ . 



Solving Eqs. (23), (24), (26), (31) and (32) numerically, we obtain the time evolu- 



tion of scalar fields, the energy densities, the scale factor and the lepton asymmetry. 



_ \2 /\-f^ X^ Ay"2 \2 ?\,y2 ■, i , 

Note that limjv^o r/f-s-AfiWi — sItM'^ (-^ ^ "^ k?m'^ )(-^ ^ '^ iSm'^ ) ' ■ The expressions in brackets 
emphasize that the decay is kinematicaUy possible if m\j' > 2m)^'. However in the parameter range 



of interest, Xu/Mp <C n/M holds, as can be seen from Eqs. (20). Thus we have Tn^NiNi — g^rA/a • 

*The assumption T^i < ^n^ (see section^ implies A33(j/i/yJ)33 > Aii(y^yJ)ii. The assumption of 
hierarchical heavy neutrinos implies A33 ^ An. 
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respectively. The former two are displayed in Fig. [3j The regime of reheating is char- 
acterized by oscillating scalar fields and can be divided into distinct phases: After the 
end of inflation both A^ and H fall to their true minimum and begin to oscillate. After 
only a few oscillations the classical field H settles at its minimum and the dynamics of 
the system is governed by the oscillation of the A^ field. The further evolution of the A^ 
oscillations is governed by Hubble damping. As long as "H ^ V^ the universe is gov- 
erned by (damped) oscillating scalar fields which can be interpreted as (decaying) heavy 
particles. This implies a matter dominated universe out of thermal equilibrium. Ultra- 
relativistic particles are produced through the decays of the heavy particles, however 
they are diluted by the expansion of the universe. As soon as 7{ ~ F^r the radiation 
energy density becomes dominant and the light particles begin to thermalize. This 
marks the end of reheating and determines the reheat temperature and the asymmetry 
ul/s. 

4.2 Simplified Treatment with Boltzmann Equations 



Since the set of equations (23), (24), (26) and (31) is quite involved, a common attempt 



in the literature (e.g. [36]) is to simplify these equations by time-averaging the equations 
of motion of the scalar fields. The result is a set of Boltzmann equations for the matter 
energy density p^ = Pn + Ph and the radiation energy density pr completed by the 
Priedmann equation 



3 Ml 



Pm + ^'HpM — —'^nPm, 
P'r + 4:1-1 PR = TnPm, 

{Pm + Pr) = T-i^- 



(33) 
(34) 

(35) 



The lepton asymmetry is determined by Eq. (32). The big advantage is that these 



equations have approximate analytical solutions. However their derivation (see e.g. [36] ) 
implies an important assumption concerning the scalar potential V{N, H). In order to 
rewrite the time-averaged kinetic energy density in terms of the total energy density by 
exploiting the Virial theorem we must assume that we can write the scalar potential as 
V{N, H) = Vn{N) + Vh{H) with V"^ ~ A^'^ and Vh ~ H"-. Eqs. (|33| - (|35|) are obtained 
with r = 2. Numerical simulations of the full system (23), (24), (26), (31) and (32) 



show that this assumption is not justified in the early oscillation phase in the model 
described by Eq. ^ since the large oscillations of the A^-field result in a highly coupled 
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Figure 3: (a) Field dynamics of the sneutrino field (blue) and the waterfall field (black), (b) Evolution 
of the vacuum (green), matter (blue) and radiation (red) energy densities. The parameters chosen for 
these plots are k = 0.5, M = 0.0032 Mp, ra^ = 4.9 • lO^^ GeV, yi = 10-^ ksh = 0.5. The timescale 
is given in units of the Planck time, tp — h/Mp ~ 2.7 • 10~^'^s. 

system with higher orders terms in the scalar potential playing a non-negligible role. 
However they do also show that for t ^ T^^ ^ "H^^ the simpler system of differential 



equations (33) - (35) does give a good approximation. This is the point of time relevant 



for the predictions of the reheating phase. 

Having seen that the results of the numerical solutions to the full field equations 
for t ^ r~^ can be approximated reasonably well by the simpler set of Boltzmann 
differential equations, we can now find approximate analytical solutions to the latter 
and use these expressions to find estimates for the reheat temperature and the produced 
baryon asymmetry 



T. 



BH 



9 



An^g, 



1/4 



y {yuyt)nmNMp , 






(to) ~ 3.45 



C 



C-1 



9;'^'e 



{yuvhii 

m]sf/Mp 



(36) 



(37) 



Combining (36) and (37) reproduces the familiar relation ub/s 



~ e 



TRH/niN (see 

e.g. |31 ES])- These results must be compared with existing bounds on the reheating 
process. The WMAP 7 year data combined with measurements of the baryon acoustic 
oscillations and todays Hubble parameter imply ^ = (6.19 ± 0.15) ■ 10~^° |37], thus 
yielding Trh > 1-4 ■ 10^ GeV. 

Additionally, the reheat temperature is bounded from above by the so-called grav- 
itino problem [Ml ESI SDl SIl 112] • A high reheat temperature would result in an over- 
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production of gravitinos. If these are stable, then the fact that their energy density can 
not be larger than the present total energy density of the universe leads to a bound 
on the reheat temperature in terms of the gravitino mass WI3/2. On the other hand, 
if gravitinos are not stable, they can either decay before or during and after the Big 
Bang Nucleosynthesis (BBN). In the former case (i.e. heavy gravitinos), with R-parity 
conserved the gravitinos will decay into the lightest supersymmetric particle (LSP) and 
their production is thus constrained by the dark matter abundance. This yields a fairly 
model independent bound of T^h < 2 ■ 10^^ GeV for an LSP mass of about 100 GeV to 
150 GeV. In the latter case (i.e. light gravitinos), the decay of the gravitinos would alter 
the outcome of BBN and create a conflict between BBN predictions and observations. 
This yields even stronger, however model dependent, constraints on the reheat tem- 
perature. Combining these arguments yields a constraint on the reheat temperature of 
typically Trh < 10^ — 10^° GeV, depending mainly on the model under consideration 
and on the value of 1713/2- The resulting preferred region in {niN, yi) -parameter space 
is depicted in blue in Fig. |4} 

4.3 Remarks on Preheating 

Note that throughout this chapter we have focussed on the evolution of the homoge- 
neous fields A^ and H. It has been pointed out that under certain circumstances this 
might not be sufficient, since modes with /c 7^ of all fields in the model can be strongly 
excited at the end of infiation and before the beginning of reheating in a process re- 
ferred to as preheating. There are two types of preheating worth mentioning in the 
context of hybrid inflation, namely preheating via parametric resonance [301 SSI IIH SS] 
and tachyonic preheating [l6l S?]. In the former case, the coupling of fermions and 
bosons to the oscillating inflaton field results in oscillating mass terms for these parti- 
cles. Solving the respective equations of motion (roughly the equation of an harmonic 
oscillator with an oscillating mass as described, e.g., by the Mathieu equation) can yield 
explosive particle production. However, in the region of parameter space of interest to 
us, any heavy particles that are produced by this mechanism will decay back into heavy 
(s) neutrinos or into radiation. The radiation produced directly or indirectly through 
this process at the beginning of the reheating phase will however be strongly diluted 
during the ongoing matter dominated phase and thus be insignificant for the outcome 
of the reheating phase. Thus in our model, parametric resonance will not affect the 
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results discussed above, mainly due to the structure of the mass spectrum and the very- 
small effective Yukawa coupling iji. 

Tachyonic preheating occurs when the squared mass of the H field becomes negative, 

(s) 
triggering the waterfall ending inflation. Modes of the H field with k < \mH \ grow 



exponentially q causing particle production of bosonic and fermionic fields coupled to 
the waterfall field [46] and creating an inhomogeneous field H{x, t) which can cause 
the formation of topological defects when the waterfall occurs [17]. It was stated 
in |17j that the production of fermions and bosons coupling to the waterfall field with 
a coupling strength g is suppressed by Pb,f/ Pv ~ 10^'^ g with pv denoting the energy 
density during infiation. Thus in the parameter region of interest, this is negligible 
in our model. On the other hand, the production of topological defects could indeed 
dominate the evolution of the universe in an early stage. However, since we have not 
observed any topological defects yet, a mechanism to prevent or dilute these objects 
(e.g. a preferred waterfall direction or a slight shift of the potential energy of the 
discrete vacua) is typically implemented. We will assume that the higher dimensional 
operators denoted by dots in Eq. ([l]) provide such a solution so that at some time after 
the waterfall, the universe is dominated by the lightest right-handed sneutrino. The 
evolution from this point on is correctly described by the classical theory of reheating, 
as discussed above. Other possible scenarios in which the evolution of the universe may 
not be dominated by the homogeneous component of the infiaton field remain to be 
explored in this context. 

5 Summary and Conclusions: Combining Inflation 
and Leptogenesis 

In sections |3] and |4] we have investigated the conditions under which infiation, with 
primordial perturbations in accordance with the latest WMAP results, as well as suc- 
cessful leptogenesis can be realized simultaneously in simple models of sneutrino hybrid 
infiation as outlined in section [2] The combined results are summarized in Fig. |4j 



^It was pointed out in jJHj that in some hybrid infiation models a fragmentation of the infiaton 
condensate can occur, causing the evolution of the universe to be dominated by these 'lumps' instead 
of by the homogeneous component of the infiaton field. However this 'lump' formation requires a 
flatter than cfp' potential (with = {N^ H}), which does not appear in our model as can easily be 
checked from Eq. (|8|. 
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Figure 4: Preferred region (95 % CL) in (totv, yi)-paranieter space from inflation and leptogenesis. The 
preferred region derived from inflation, in particular from the WMAP constraints on the spectral index 
Us, is marked in grey. The favored region obtained from reheating and leptogenesis is depicted in blue. 
The lower bound corresponds to the baryon asymmetry measured by WMAP assuming a maximal 
CP-violation ei. The gravitino problem imposes a (model dependent) upper bound on the reheat 
temperature Trh- This yields the different shadings, corresponding to a different reheat temperature 
and correspondingly a different CP-violation. Finally, the red lines depict a constant effective neutrino 



mass parameter rhi. The upper line corresponds to a light neutrino mass of 0{^/Arn^^^ ^^j), whereas 
the dashed line depicts the borderline between predominantly thermal and nonthermal leptogenesis. 
For these plots we chose k = I. (a) shows the situation for ksh — 1 = — 1. Note that after reaching 
a minimum at Ug — 0.98, the spectral index acquires large values for increasing tun. (b) depicts the 
situation for Kg^ ^ 1 = +1- In this case, the sign of the SUGRA corrections flips and Ug decreases for 
large values of rrijv • 
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The dynamics of inflation is governed by the scale M of the phase transition ending 
hybrid inflation, the mass of the lightest right-handed (s) neutrino mTv, the vacuum 
energy parameter k (= waterfall field self coupling) and the parameter k^sh control- 
ling the SUGRA corrections. In principle, terms depending on the neutrino Yukawa 
coupling matrix could contribute, too. However, in our case the comparatively small 
first generation Yukawa couplings make these contributions negligible. With M fixed 
by the amplitude of the scalar CMB fluctuations and n ~ 0(1), the spectrum of the 
CMB fluctuations is primarily dependent on the lightest right-handed (s)neutrino mass 
niN- For large values of niN, SUGRA corrections controlled by ksh become important, 
with the sign of these contributions depending on the sign of k,sh ~ 1- For the spectral 
index, its running and the tensor-to-scalar ratio the predictions are shown in Fig. [2] 
Recent WMAP observations constrain the preferred region for the spectral index n^, 
thus imposing a constraint on the preferred region for ttin- For k = 1 and Ksh — 1 = ±1 
this is marked in grey in Fig. |4j 

On the other hand, the decisive quantities of reheating and leptogenesis, namely 
the reheat temperature Trh and the baryon asymmetry UB/n^ depend on the effective 
first generation neutrino Yukawa coupling yi, the CP asymmetry ei and the mass of the 



lightest right-handed (s)neutrino rri]^ (see Eqs. (36) and (37)). The latter parameter 
is thus the link between inflation and leptogenesis. The preferred region of parameter 
space resulting from bounds on these quantities is marked in blue in Fig. |4} It is 
bounded from below by the experimental value of the baryon asymmetry measured by 



WMAP and by an upper bound on the CP-violation per (s)neutrino decay Eq. (28). 
From above it is bounded by constraints imposed on the reheat temperature from 
the gravitino problem. Since these are model dependent, we have plotted the regions 
satisfying Trh < 10^°, 10^, 10^, 10^ GeV in different shadings. Note that a higher 
reheat temperature at a fixed value for rriN automatically corresponds to a smaller 
value of ei in order to match the measured baryon asymmetry. The resulting preferred 
region in parameter space implies an effective first generation Yukawa coupling yi = 
(9(10^^ — 10"'^). The upper part of this range is of the same order as the first family 
quark and charged lepton Yukawa couplings, which in the MSSM with moderate tan/3 
are of the order 10~^ to 10~^. 

Throughout this paper, we have assumed nonthermal leptogenesis and hierarchi- 
cal masses of left-handed as well as right-handed neutrinos. Assuming that the light 
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neutrinos obtain masses via a type I seesaw mechanismQJ both assumptions depend 
on the value of the effective hght neutrino mass parameter (also dubbed washout pa- 



rameter) rhi = y'liv)'^ /ui]^. More explicitly, one can easily see from Eq. (36) that 
{Trh/^^nY ~ 4.0 ■ lO^mi/eV. Lines of constant rhi are marked in red in Fig. El cor- 
responding to a fixed relation between Trh and mAr. Simultaneously, they give the 
order of magnitude for the mass of the left-handed neutrino rriy^ ~ rhi. In the pre- 
ferred region of parameter space, we find rhi < 3.4 ■ 10~^ eV, thus implying nonthermal 



leptogenesis with Trh <^ mN and ruui ^ \/ ^"^atm, soi- 

Finally, we want to comment on possible extensions of this scenario and the signif- 
icance of cosmological observations in the near future. In Fig. |4] we have set k = 1. 
Allowing for 0.5 < k < 2 gives qualitatively the same picture (see section Isl in par- 
ticular Fig. PI with a somewhat shifted grey region. For example, for k = 2 the grey 
region is extended to the left to rriN = 2 ■ 10^*^ GeV whereas for k, = 0.5 it is extended 
to the right to m^ = 7 ■ 10^^ GeV for ksh ~ 1 = +1o Another interesting possibility 
would arise if the experimentally preferred region for the spectral index was raised, fa- 
voring a spectral index closer to 1. This would lower the preferred range for the lightest 
(s) neutrino mass rriN significantly and thus open up the region of thermal leptogenesis 



and allow for m,^-^ ~ 0{.^/Aml^^^ ^^j). The forthcoming results of the Planck satellite 
will make the requirements for rrii^ more accurate. 

In summary, we have pointed out that successful sneutrino hybrid inflation and 
leptogenesis can be achieved in this framework, and that combining both imposes re- 
quirements on the parameters of the underlying particle physics model. We obtain a 
mass for the lightest right-handed (s)neutrino of mj\f = (9(10^° — 10^^) GeV, an effec- 
tive first generation neutrino Yukawa coupling yi = O{10^^ — 10^^) and a very light 
left-handed neutrino with m,y-^ < O{10~^) eV. Furthermore, we find that leptogene- 
sis occurs via nonthermal leptogenesis (with Trh/^^n < 0.1 for /€ = !). Concerning 
the spectrum of the CMB fluctuations, we predict a running of the spectral index of 
—0.0004 < as < 0.0002 and a tensor-to-scalar ratio of r < 0.015. Our results pro- 
vide a guideline for the construction of explicit particle physics models incorporating 
sneutrino hybrid inflation and subsequent baryogenesis via nonthermal leptogenesis. 



^''This implies a mass matrix for the left-handed neutrinos {my)ij = ^{y^M ^yi,)ij{v)'^ /2. 
^^For Ksh = —1 the respective region is extended to niM = 2 • 10^^ GeV. 
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